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SUMMARY 


One year of longwave radiation data from July 1975 through June 1976 from 
the Nimbus 6 satellite Earth radiation budget experiment is analyzed by repre- 
senting the longwave radiation field by a spherical harmonic expansion. The 
data are from the wide-f ield-of-view instrument. Results show that the limit 
of the spherical harmonic representation is 12th degree, based on degree 
variance plots from 12 months. Degree variance plots also show that most of 
the power is in the lower degree terms. The axisymmetric (zonal) terms domi- 
nate, with their coefficients representing approximately 80 percent of the 
degree variance. Contour maps of the radiation field for 12 months show the 
geographical distribution of Earth emitted radiant exitance (W-m"^) and reveal 
areas of high and low emitted radiation. The analysis also shows differences 
in the Earth emitted radiation between the Northern and Southern Hemispheres, 
presumably due to land and ocean distributional differences. 


INTRODUCTION 

The nonuniform distribution of radiative heating and cooling over the 
Earth drives the motion of the atmosphere and oceans, and an understanding of 
the interaction between the radiation and the dynamics of the atmosphere and 
oceans is required for understanding climate processes. For this reason, 
several Earth radiation instruments of increasing sophistication have flown on 
spacecraft. As early as 1960, instruments for measuring the Earth’s radiation 
budget were flown on Explorer 7 (Vender Haar (1968)). Since then, a number of 
satellites have made Earth radiation budget measurements. The ESSA meteoro- 
logical satellites, the Nimbus satellites, NOAA 1 , and ITOS 1 are some of the 
satellites which have made Earth radiation budget measurements. Radiometers 
flown on these satellites have provided a wealth of Earth radiation budget 
data, including longwave measurements, shortwave measurements, wide-f ield-of- 
view (WFOV) measurements, and scanning radiometer measurements. Descriptions 
of these satellites, and analysis of data from them, have been reported by 
numerous authors such as Raschke et al- (1973), Weaver and House (1979), 

Ellis et al. (1978), and Green and Smith (1978). 

The most recent instrument flown on a satellite for making Earth radiation 
measurements was the Earth radiation budget (ERB) instrument aboard the Nimbus 6 
and Nimbus 7 spacecraft (Smith et al. (1977); Jacobowitz et al. (1979))- The 
ERB instrument includes two Earth viewing WFOV radiometers, one which measures 
total ir radiance (0.2 to 50 \im) and another which measures shortwave (SW) 
irradiance (0.2 to 3.8 ym) . The difference between the two irradiances is the 
longwave (LW) or Earth emitted radiation. The ERB instrument has been pro- 
viding WFOV measurements since July 1975. These WFOV data are suitable for 
studying large-scale processes over large spatial scales. WFOV sensors inte- 
grate radiant exitance (W-m“^) from all portions of the Earth atmosphere systems 
within their fields of view. For the Nimbus 6 spacecraft at an altitude of 
about 1100 km, the field of view is the entire disk of the Earth. 



A deconvolution technique (Smith and Green (1976); Green and Smith (1978)) 
is applied in this paper to enhance the resolution of the WPOV measurements. 
Deconvolution represents the radiant exitance at the top of the atmosphere 
(denoted by TOA and defined herein to be at an altitude of 30 km) by a trun- 
cated series of spherical harmonics. Much of the earlier work on radiation 
budget from satellites has used the geometric shape factor technique to give 
radiant exitance at TOA. The two techniques are equivalent for global average 
estimates, but deconvolution gives more information at higher orders. 

The Nimbus 6 ERB WFOV data have been analyzed by various investigators. 
Jacobowitz et al. (1979) recently analyzed and published the results of the 
first 18 months of ERB measurements. The inverse square approximation (geo- 
metric shape factor) was used to determine the radiant exitance at TOA. Some 
of the earlier investigations have concentrated on the data for August 1975. 

One such analysis, reported by Smith et al. (1977), defined the monthly mean 
radiation field by a global contour plot and the global mean. The data anal- 
ysis technique also employed the inverse square approximation. Green and Smith 
(1978) analyzed the same data using the technique of parameter estimation and 
the concept of deconvolution. The parameter estimation technique used in the 
analysis by Green and Smith is based on solving a system of simultaneous 
measurement equations using least squares. 

In the present paper, a deconvolution (resolution enhancement) technique 
is presented to study the annual cycle of the distribution of the Earth's 
LW radiation. The study is based on 1 year of WFOV data from July 1975 
through June 1976. The resolution enhancement technique was used to analyze 
the data and extract mean radiation fields on a monthly basis. However, 
instead of considering each measurement individually, the data were averaged 
over 5^ by 5^ grid areas at satellite altitude. This averaging technique 
reduced the computational burden at the expense of smoothing the data. This 
data set was analyzed to produce radiant exitances at TOA. The deconvolution 
technique takes advantage of the fact that spherical harmonics are the eigen- 
functions of the measurement operator and reduces the radiant exitance field 
from satelite altitude to TOA by dividing by the appropriate eigenvalues. 


SYMBOLS 

aJJ, bJJ complex coefficients of spherical harmonics for degree n and 
order m, W-m“^ 

Cj^ integral defined by equation (14) 


^n' ^n real coefficients of spherical harmonics at top of atmosphere 
for degree n and order m, W-m“2 

'^m '^m 

Cj^, Sj^ real coefficients of spherical harmonics at satellite altitude 
for degree n and order m, W-m*“^ 


da 


differential surface element 
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angular response function 

satellite altitude above r^, 1070 km 

integral defined by equation (16) 

total number of global 5^ by 5^ grids, 1654 

grid region 

reflection about equator of kth grid region 
r adi ance , W-m” sr “ ^ 

linear integral measurement operator defined in equation (4) 
radiant exitance at top of atmosphere, W-m”^ 
measurement, W-m“^ 

average of measurements in kth grid region, W-m"^ 
degree of spherical harmonic expansion 

normalizing factor for spherical harmonics of degree n and 
order m 

associated Legendre polynomials of degree n and order m 

limb-darkening function for emitted radiation 

radius to top of Earth atmosphere system, 6408.165 km 

integral defined by equation (15) 

time 

complex spherical harmonics of degree n and order m 

real spherical harmonics of degree n and order m (see 
eqs. (9) and (10)) 

cone angle at satellite from nadir to surface element at top of 
atmosphere, deg 

Earth central angle between satellite nadir to a point on surface 
at top of atmosphere, deg 

Kronecker delta 

colatitude, deg 
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©g cx5latitude of subsatellite point, deg 

6 zenith angle of exiting ray, deg 

nth eigenvalue of measurement operator 
degree variance for degree n, 

$ longitude, deg 

longitude of subsatellite point, deg 
({) azimuth angle of exiting ray, deg 

Q solid angle at satellite subtended by surface element, sr 


DATA SAMPLE 

The Earth radiation budget (ERB) instrument aboard the Nimbus 6 satellite 
obtained Earth radiation measurements with wide-angle and scanning narrow-angle 
radiometers- A description of the ERB instrument is given by Smith et al. 
(1977)- Only measurements with the wide-angle radiometers are considered in 
this paper- The wide-angle radiometer sensed radiation in two channels; one 
measured total irradiance (0.2 to 50 ym) , the other, shortwave irradiance 
(0.2 to 3.8 ym) . Twelve months of these measurements, beginning July 2, 1975, 
and ending June 29, 1976, were supplied on magnetic tapes by the National 
Oceanic and Atmospheric Administration (NOAA) . Only the longwave radiation 
component is analyzed herein, and this component was obtained by subtracting 
the measured shortwave component from the measured value of total irradiance. 
All measurements from the data tapes are at satellite altitude; and in the 
analysis which follows, all measurements are increased by 11 percent as a 
calibration correction. This correction is discussed by Smith et al. (1977) 
and accounts for a discrepancy between the ERB experiment fixed wide-angle 
measurements and the integration of scanning narrow-angle measurements. 

For most of the data period (July 2, 1975, to June 29, 1976) the ERB 
instrument operated on a nominal duty cycle of 2 days on and 2 days off; how- 
ever, because of operational considerations, this schedule was not rigorously 
followed. Measurements were taken at 4-second intervals during the time the 
instrument was operating. However, the data tapes supplied by NOAA and on 
which the analysis in this paper is based reduced the amount of data by 
averaging every four measurements. This caused some smoothing and gave 
irradiance measurements at 16-second intervals. Figure 1 is a 12-month 
calendar showing the approximate time periods covered by the data set as 
received from NOAA. Each day on the calendar begins at 12 midnight, and the 
dark bands represent the approximate time periods during which data were 
recorded. 
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DATA EDITING 


Some measurements were edited from the data set because near sunrise and 
sunset the WFOV sensor was exposed to the Sun (Jacobowitz et al. (1979)). 
Exposure to the Sun was due to the design of the wide-angle radiometer. In 
order for the wide-angle radiometer to view the entire Earth disk, the actual 
field of view is about 8^ larger than the Earth disk and includes a thin annulus 
of deep space. When the satellite is near sunrise or sunset, the Sun comes into 
view and contaminates the measurements. When the Sun is on the horizon at sun- 
rise and sunset, the Sun zenith angle at the subsatellite point is 121.5^. This 
is shown in figure 2. Allowing for a 2^ margin on each side of the thin annulus 
resulted in a range of Sun zenith angles from 111.5^ to 123.5^. Based on plots 
of SW and LW measurements near sunrise and sunset for a few orbits during the 
month of August 1975, deleting measurements when the Sun zenith angle ranges 
from 111.5*^ to 123.5^ is sufficient for eliminating Sun contaminated 
measurements . 

The Sun contaminated measurements introduce an additional sampling bias in 
the Nimbus 6 data. (A bias in the data already exists because the orbit of 
Nimbus 6 is Sun synchronous, meaning that the radiant exitance at any lati- 
tude is based on measurements made only at two local times.) Since the Sun 
contaminated measurements occur at night (Sun zenith > 90^) , the editing out 
of these data reduces the number of samples at night while leaving the daytime 
data set intact. If the day and night measurements are significantly different, 
then the radiant exitance for the zone in question will be biased toward the 
daytime values. No attempt was made to adjust for this effect. 

The zone showing the deficiency in number of measurements depends on the 
time of year. Figure 3 shows the number of measurements for different zones 
for the month of September. The maximum deficiency occurs between ±(60^ to 
65^) latitude. Figure 3 also shows that the least number of measurements 
occurs in the ±(80^ to 85^) zone. The measurements in this zone were all 
located near ±80*^ since the orbit only goes to about ±80.1^ latitude. The 
largest number of measurements occurs in the ±(75^ to 80^) zone, and the zones 
in between are sampled about equally except for the deficiency due to Sun 
contaminated measurements. The number of zonal measurements is dictated by 
the inclination of the Sun synchronous orbit of Nimbus 6. The dashed line in 
figure 3 shows how the number of measurements would change if Sun contaminated 
measurements were not deleted. Between 6 and 7 percent of the measurements for 
September were deleted because they were contaminated by the Sun. It should 
also be noted from figure 3 that when Sun contaminated measurements are not 
deleted, the number of measurements at about 65^ north latitude is less than 
at 65^ south latitude. This data dropout is due to the data acquisition sta- 
tion in Alaska. 

A possible source of error in the longwave measurements which could vary 
from 1 to 4 W-m“^ was caused by a thermal transient in the shortwave measure- 
ments due to irradiation of the sensor at spacecraft sunset (Green and Smith 
(1978)). No attempt was made to correct for this transient in the data. 
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II 


I I I 


Measurements were also omitted if they were greater than 240 W-m”2 or less 
than 50 W-m"*^ at satellite altitude. The criterion for keeping measurements 
between the range 50 and 240 W-m“^ was based upon visual inspection of plots 
showing radiant exitance versus time for all measurements. Another criterion 
deleted a measurement if its value changed by more than 10 W~m**^ from the pre- 
ceding measurement over a 16-second interval. This criterion is based on the 
decision to omit measurements with values greater than the 3a value of the 
difference of all adjacent measurements, where adjacent measurements are not 
separated in time by more than 16 seconds. The 3a value for the 3-month period 
of July, August, and September is 8 W-m' 


-2 


Parts of the 12-month data set had anomalous measurements even after the 
preceding editing had been completed. The anomalies were discovered by visu- 
ally inspecting radiant exitance-time plots of all measurements for each month 
of the 12-month data period. Figure 4 shows a composite of one of these plots 
which shows some anomalies in the data. The anomalies are characterized by 
periods of unusually low, unusually high, and erratic radiant exitance values. 


Measurements for a 2-hour period on July 10 and a 2.5-hour period on 
July 17 were deleted because of erratic spiked data. About 2 days of measure- 
ments were deleted between July 21 and July 23 because average values suddenly 
dropped to about 90 W-m”^. The reason for this sudden decrease in radiant exit- 
ance is not known. Another questionable set of data was omitted on August 16, 
1975. Starting at 1000 hours GMT, the measurements gradually decreased in value 
during a full satellite revolution and then stabilized at values approximately 
30 W-m“^ lower than expected. This problem occurred at the end of a data acqui- 
sition period; however, the measurements from the next data period did not 
exhibit the problem. The reason for the gradual decrease in radiant exitance 
is not known. Approximately 15 hours of measurements were deleted between 
September 5 and September 6 because of a sudden increase in radiant exitance. 
Measurements were deleted from October 23 through October 31 because of opera- 
tional problems. The net effect was incorrect radiant exitance averages for 
some grid regions. About 3 hours of data in January 1976 were edited out 
because of anomalous measurements on January 12. Sixteen hours of measurements 
on February 19 and 30 hours of measurements between February 28 and February 29 
were edited out. Measurements taken during the first half of March wec»^ edited 
out. The February and March measurements were deleted because of operational 
problems. Ten hours of measurements on June 14 were edited out because of 
erratic values of radiant exitance. 


After anomalous measurements were deleted from the data, a few question- 
able points still remained in the data. The following procedure was used in 
an effort to delete as many of these anomalous measurements as possible. 
Measurements were grouped into 5° latitudinal bands. Any measurement within 
a band that differed by more than 2a from the average radiant exitance in that 
band was edited out. This editing procedure caused about 4 percent of the 
measurements to be deleted. 


The total number of measurements for the 12-month period, on which the 
analysis in this paper is based, was about 588 000 after all editing. This 
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figure represents about 30 percent of the maximum number of m^iasurements for 
the year had th^ ERB instrument made measurements at 1 6-second intervals 
uninterrupted. 


ANALYSIS TEX3HNIQUE 

The ERS wide-angle data were used to estimate the longwave radiative field 
at the top of the atmosphere (defined herein as 30 km) . The radiative field 
was represented by a spherical harmonic expansion, and the coefficients were 
estimated with measurement data. 


Measurement Model 

The top of the Earth atmosphere system was approximated by a sphere of 
radius r^. The emitted radiance L leaving any point E on this surface is, 
in general, a function of colatitude 0, longitude $, zenith angle 0, azimuth 
angle (j), and time t. (See fig. 5.) Thus, the measurement m at satellite 
altitude h, colatitude longitude Og is given by 

m(0g,^s,t) = i L(0,$,0,ct),t) g(a) d;^ (D 

^FOV 


where Q is the differential of solid angle at the satellite subtended by the 
surface element, a is the cone angle at the satellite from the local vertical 
to the surface element, and the integration is over the field of view (FOV) of 
the sensor (Green and Smith (1978)). The function g (a) is the angular 
response of the sensor to incoming radiation and, for the ERB fixed wide-angle 
radiometer, has been modeled as a perfectly black flat-plate sensor, or 


g(a) = cos a 


( 2 ) 


Also, the radiance field was modeled as being independent of (j>; that is. 


L(0,<^,0,t) 


1 

7T 


M(0,$,t) 


R(e) 


(3) 


where M is the radiant exitance at the top of the atmosphere and R(0) is a 
directional function assumed to be independent of position and dependent only 
upon the zenith angle of the exiting ray. The time dependence is removed by 
considering monthly averages of measurements. Substituting equations (2) 
and (3) into equation (1 ) yields 
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m(©s,$s) 


M(©,$) R(9) cos ot dfl 


(4) 



or 


m(0s/'ts) =^[m(0,$)] 


where if denotes the linear integral measurement operator of equation (4). It 
has been shown by Smith and Green (1976) that the eigenfunctions of this linear 
operator are spherical harmonics for any R(9); that is. 


^ [^nj (0s''^s^ “ ^n ^n(®s»‘^s) (5) 

where Yn(©s»'^s) is a spherical harmonic of order m and degree n evaluated 
at the subsatellite point. The associated eigenvalue is given by 


^n = 


'^a=o 


h 0, 

Pn(cos 


Y) R(6) cos a sin a da 


( 6 ) 


where Pn(cos Y) denotes the Legendre polynomial of degree n as a function 
of the Earth central angle Y» (See fig. 5.) The integration is from nadir 
(a = 0) to aj^, where aj^ is the cone angle to the horizon. Eigenvalues 
through degree 12 are listed in table I for both a Lambertian and limb-darkening 
model (Raschke et al. (1 973); Green and Smith (1 978)). 

Over most of the region subtended by the Nimbus 6 WFOV radiometer, g(a) 
has a cos a response except for an enhanced response near 60® which is thought 
to be due to multiple reflections from the ERB instrument. However, since the 
amount of energy contributed to the measurements from 60° is small, it is felt 
that modeling the radiometer as a cosine response is adequate. 

The effect of choosing one directional function for the entire globe has 
been investigated. An extreme case was considered where the entire globe was 
assumed Lambertian, then the entire globe was assumed greatly limb-darkened and 
the two results were compared. The results reflect the sensitivity of R(0) 
being independent of position. The two solutions gave identical global means. 
The area-weighted mean of the absolute zonal difference of the Lambertian solu- 
tion and the limb-darkened solution was 0.55 W-m“^. The area-weighted mean of 
the absolute 1 0° regional difference of the Lambertian solution and the limb- 
darkened solution was 0.75 W-m”2. The directional function used in the present 
study (Green and Smith (1978)) is intermediate to the two directional functions 
discussed above and is considered to yield good results. 
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Deconvolution 


Because spherical harmonics are eigenfunctions of the measurement operator, 
it is convenient to represent the radiant exitance at TOA by a truncated series 
of spherical harmonics; that is. 


N n 
n=0 m=-n 


M(0,$) 


■I 


Also, let the measurements be represented as 


N n 


m(©,<D) = ^ Yn(0, 




n=0 m=-n 


and from equations (4) and (5) it follows that 


N n m 
ar 


M(0,<i)) = Y, Y. r 




n=0 m=-n 


W-m 2, are complex coefficients of spherical harmonics. 


Estimation of Coefficients 

For the purposes of this analysis, a real formulation of the spherical 
harmonics was used. Thus, the radiant exitance at the top of the atmosphere 
is 


N n 

M(0,$) = Y Y + Sn Yn(0,$)] (7) 

n=0 m=0 


where 


m 


and bi 


m 

h f 


in 
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''in ^ni 

where and are real coefficients and are estimated from the measure- 

ments such that 


n 


m(0,$) = ^ ^ Yn(0,$) + S™ yJ(0,$) 


(8) 


n=0 m=0 


Moreover, the following definitions apply: 


Y™(0,$) = Nn cos(m$) Pn(cos 0) 


(9) 


yJJ(0,$) = nJJ sin(m$) P™(cos 0) 


( 10 ) 


n"* = 


(2n + 1 ) (n - m) ! 


(n + m) ! 


1/2 


01 ) 


X 


vm 




da 


Sphere 


= I 


-m 


'vi; 


h . k 


da = 4tt 6n 5 


n 


Sphere 


( 12 ) 


X yS 5^ ® . 0 

Sphere 


( 13 ) 


Also, the following notation will be convenient: 


Cm(^l»'^2) ) cos(m$) d$ (14) 
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Sm<^l »'^2) 


(15) 




sin(mO) 


d$ 


p 0=©2 

(©1 ,© 2 ) = \ P^(cos 0) sin 0 d0 

"^0=01 


(16) 


Multiplying equation (8) by y 5^(0,<|)) and integrating over the sphere, taking 
into account the orthogonal conditions of equations (12) and (13), gives the 
coefficients at satellite altitude as a function of the measurements; that is. 



{4ir)-l \ m(0,$) Yn(0,$) da 

^Sphere 


(17) 


However, the measurements m(0,<^) are not known everywhere but only at dis- 
crete points, so that equation (17) cannot be applied directly. This problem 
was overcome by dividing the sphere into grid regions and assuming m(0,$) 
was constant over each region. Thus, mj^ denotes the average of all measure- 
ments that are in the kth region. The coefficients in equation (17) are given 
by 


C 


m 

n 



mk 


I 


TT ! 

Yn(0,$) 


Region k 


da 


and from equation (9) 


K 


*< l >==$2 W 


c?: = n^4tt)-i y m, r 

(k) 


p0=©2(k) 

cos (m^) d$ \ Py^(cos 0) sin 0 d0 

"^0=©1 (k) 


and from equations (14) and (16) 
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( 18 ) 


K 

C™ = <(47r)-‘' ^ nik C„(k) l™(k) 
k=l 


Similarly, 


K 

Sn = Nn(4TT)-‘' ^ ni,^ Sm(k) I™(k) (19) 

k=l 


The regions as defined in this analysis are approximately equal area 
regions, the sides of which are latitude and meridian lines. At the equator 
the regions are square. For 5^ by 5^ regions, the sphere is divided into 
5^ latitudinal bands and each band is divided into an integral number of areas 
such that each area has approximately the same area. The band from 0^^ to 5^ 
colatitude contains three regions. The band from 5^ to 10^ colatitude con- 
tains nine regions, etc., down to the band from 85^ to 90^^ colatitude, which 
contains 72 square regions, 5^ by 5^. The total number of regions for a 
50 50 gj:i (3 system is K = 1654 for the entire globe. 

Equations (18) and (19) provide a very efficient way to determine the 
spherical harmonic coefficients of the measurement representation. First, the 
total number of measurements is greatly reduced by averaging all measurements 
in the same region. The integrals Cj^ and are rapidly computed recur- 

sively with the trigonometric addition formulas. Recursive formulas for 
are given in the appendix. The computational effort is halved by taking 
advantage of the symmetry of the grid system about the equator. If k* corre- 
sponds to the region that is the reflection about the equator of the kth region, 

then it can be shown that lJJ(k*) = iJJ(k) ./ Thus, k in equations (18) 

and (19) need only range over the regions in the Northern Hemisphere, and 

K/2 

Cn = Nn(47r)-1 ^ C„,(k) In(k) + Cm(k') In(k*)J (20) 

k=l 


or 
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K/2 


Cn = N?J(47r)-'> ^ [mk + (-1 ) C„(k) I™(k) 

k=l 


Similar ly. 


K/2 


im 


= Nn(4Tr)-‘> ^ [n>k + (-1)^+1%,^, j Sm(k) In(k) 

k=l 


(21) 


The spherical harmonic coefficients of the measurement representation at 
satellite altitude are computed using equations (20) and (21). The coefficients 
for the top of the atmosphere are given by 


C 


m 

n 


i.m 


and 


( 22 ) 


S 


m 

n 


^m 

Sn 


(23) 


where is the nth eigenvalue from equation (6) based on the limb-dar kening 

function R(6). When referring to the time-dependent nature of these coeffi- 
cients, C^(t) and S^(t) are used. 


RESULTS AND DISCUSSION 

One year of longwave radiation data from the ERB instrument aboard Nimbus 6 
and taken during the period July 1975 through June 1976 has been analyzed to 
estimate the longwave radiant exitance at the top of the atmosphere (TOA) . The 
radiative field was represented by a spherical harmonic expansion, and the coef- 
ficients were estimated with measurement data. 
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Degree Variance 


The spherical harmonic coefficients competed using equations (22) and (23) 
can be used to define the spatial spectrum of the radiation field in terms of 
degree variance. From Green and Smith (1978), degree variance for a given 
degree n at satellite altitude is defined to be 



m=0 


At the top of the atmosphere (TOA) , the degree variance for a given degree n 
is 


TOA 



The degree variance for all n is analogous to the amplitude spectrum in 
ordinary harmonic analysis and gives an indication of the power in each spa- 
tial frequency. Degree variance plots for each month as a function of degree n 
at satellite altitude and at TOA are used to establish the limit in degree to 
which a spherical harmonic representation of the data can be realized. Figure 6 
shows the degree variance plot for the average of the 1 2- month data set for each 
degree 1 through 30. The lower curve is the degree variance at satellite alti- 
tude, and the upper curve represents the degree variance at TOA. Divergence of 
the two curves begins at about degree 12. 

Figure 7, similar to figure 6, shows the standard deviation of average 
degree variance over 12 months at the top of the atmosphere for degree 1 to 

2 

degree 30. During the year there are 12 values of for each degree n. 

Standard deviation of average degree variance at satellite altitude is similar 
to the variation at the top of the atmosphere. 

The degree at which the spherical harmonic expansion should be truncated 
impacts the uncertainty of the radiant exitance contour map solutions. Based 
on a 1 2- month data set, figure 6 shows the resolution limit of the spherical 
harmonic representation to be about 12th degree, with the spectral region of 
higher frequency components beyond 12th degree produced mainly by noise com- 
prised of variations in the radiation field during the month and measurement 
errors. Monthly sampling does not appear to affect the 12th degree limit on 
truncation. This can be seen by considering the LW radiation data for 
March 1976, for which the time average of radiant exitance covered only the 
second half of the month because of some operational problems. However, the 
resolution limit was still 12th degree based on the degree variance plot for 
March. 
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Spatial sampling of the data for August 1975 averaged over the month for 

]0 -jO 

5° by 5^ grid regions and then for 2- by 2- grid regions showed that the 

2 2 

resolution limit did not change when grid size was changed. 

Attempting to represent the radiation field beyond the resolution limit 
of the sensor leads to incorrect values of radiant exitance at TOA because the 
noise is being amplified by the eigenvalues 

To give some idea of how sensitive the solution is to the degree of trun- 
cation, the measurements for the month of August 1975 were carried out to a 
15th degree expansion and then compared with a 12th degree expansion for August, 
both for a Lambertian model. Results show no difference in global mean, pole- 
to-pole gradient, and equator -to-pole gradient. The area-weighted mean of 
absolute zonal differences of 15th degree minus 12th degree was 0.85 W-m"^. 

The area-weighted mean of absolute 10^ grid regional differences was 3.94 W-m”^, 


Spherical Harmonic Coefficients 

Spherical harmonic coefficients for 12 months of longwave radiation data 
were computed using equations (22) and (23) which are based on the monthly mean 
value of the measurements in each of the 1654 grid regions. Results are for a 
degree 12 spherical harmonic expansion. For a 12th degree expansion, 169 coef- 
ficients are required to specify the radiation field. 

For a 12th degree spherical harmonic expansion approximately 80 percent of 

total degree variance is accounted for by the 12 axisymmetr ic coefficients 

to C^2* first five axisymmetr ic coefficients represent about 70 percent 

of the summed degree variance. The axisymmetr ic terms for the 12 months of 
Nimbus 6 longwave radiation data are listed in table II. Nonaxi symmetric terms 
through degree 12 and through order 12 for the 12 months are shown in table III. 

Terms below the indicated stairstep are cosine terms C^; terms above the stair- 

step are sine terms Sj^. Units of all terms are in W-m 


Geographical Distribution of Longwave Radiation 

Radiation contour maps for each month are shown in figures 8(a) to 8(1). 
Figure 9 shows monthly averages of zonal emitted radiant exitance for each of 
the 12 months. In essence, figure 9 shows a profile of each contour map where 
radiant exitance has been averaged over longitude at each 5® latitude. The 
patterns in the maps (figs. 8(a) to 8(1)) are mainly influenced by surface 
temperature, by the mean cloud cover and cloud height, and by the temperature 
lapse rate and moisture content of the atmosphere. These maps show the geo- 
graphical distribution of emitted longwave radiant exitance for the 12-month 
period and are based on a degree 12 spherical harmonic representation of the 
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radiation field. It can be seen from the contour maps, from figure 9, and from 
the table of spherical harmonic coefficients (table II) that the radiation field 
for each month is zonal in character and especially well marked in the Southern 
Hemisphere. 

The contour maps also define the variation of major features for each 
month. Centers of high and low radiant exitance show some shifting in position 
throughout the year. Highs are generally located in the tropical and subtrop- 
ical regions between -35^ latitude and +35^ latitude. A narrow band of lows in 
temperate and low latitudes between -10^ and +20^ appears throughout the year. 
These lows appear over areas with high cyclonic activity. Lows extend to high 
latitudes in both hemispheres, with the extension being a dominant feature in 
the Northern Hemisphere. Noteworthy are the lows over Greenland, Siberia, 
Antarctica, Central Africa, and Southeast Asia, the highs over Northern and 
Southern Africa, and the high over Australia. The lows over Antarctica show up 
during each month and are presumably picking up the Antarctica ice cap. Lowest 
emission throughout the year occurs over Antarctica. This low shows some shift- 
ing and variability. Emission over Antarctica is in step with the seasons, 
increasing from September through March and decreasing during the next 6 months. 
Greenland, the Beaufort Sea off the coast of Alaska, Siberia, and the Siberian 
Sea with its thick ice shield show the lowest emission over the Arctic, but not 
as low as the emission over Antarctica. In terms of temperature, the Arctic is 
much warmer on the average throughout the year than Antarctica. 

Another area of low emission is Southeast Asia during the monsoon summer, 
when the amount of longwave radiation emitted over Southeast Asia is less than 
the emitted radiation over Greenland during the same time period. This is due 
to monsoon clouds being relatively opaque with surfaces that are usually high 
and cold. During the winter months the emitted radiation over Southeast Asia 
increases because the winter monsoons bring dry, clear weather for several 
months. 

The low over Central Africa appears during each month and shows very 
little movement. This low is located in the vicinity of the Congo Basin and 
the Equatorial rain forests, which are characterized by continuously warm weather 
and abundant rainfall. During the summer season, highs exist over Northern and 
Southern Africa in the vicinity of the Sahara and Kalahari Deserts. These highs 
show some movement with the seasons. There is also a high that appears consis- 
tently near the western coast of the Australian desert region. This high shows 
some movement with the seasons. 

The monthly contour maps were shown to have a 12th order spherical harmonic 
expansion. The sensitivity caused by changing the order of expansion shows up 
clearly in the contour maps. Green and Smith (1978) analyzed the Nimbus 6 long- 
wave data for August 1975 but only produced a 5th order spherical harmonic 
expansion. Their contour map for August is qualitatively similar to the August 
contour map from the present analysis. However, the two maps differ in both 
location and spatial extent of high and low areas of radiation. In effect, the 
5th order expansion for the month of August does not contain as much high fre- 
quency information of the radiation field as does the 12th order expansion. 
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Time Histories of Zonal Coefficients 

Examination of the cJJ(t) and sJJ(t) presented in table II and table III 
shows that the axisymmetric terms (m = 0) dominate up to n = 5 with about 
70 percent of the degree variance in the first five zonal coefficients. Some 
of these coefficients have obvious physical meanings. 

Figure 10 shows the yearly cycle of the global average CQ(t) of emitted 

longwave radiation. The plot of Co(t) represented by the square symbols in 
figure 10 is the only term from the present analysis that can be readily com- 
pared with global average values of emitted radiation from Nimbus 6 published 
by other sources, such as Ellis et al. (1978) and Jacobowitz et al. (1979). In 
making the comparisons, TOA is at 30 km above the Earth's surface in all cases, 
and the radiant exitant values are averages of daytime and night ime measure- 
ments. All data sets incorporated a calibration correction referred to by Smith 
et al. (1977). Both Ellis and Jacobowitz used the inverse square approximation 
to determine radiant exitance at TOA. All three data sets show the same trend 
in global average for the year. However, the global average of radiant exitance 
based on the work by Ellis et al. (1978) is about 4 W-m"^ higher (on the aver- 
age) than the global average of radiant exitance from the present work. From 
January through June the global average from Jacobowitz et al. (1979) compares 

very well with the global averages Co(t). Also included in figure 1 0 is a 
29-month data set by Ellis et al. (1978) of global average. This data set is 
comprised of a number of sources. The trend is very irregular, with unusually 
high values of gobal radiant exitance in April and May. Some differences in 
global averages as reported by the different investigators (fig. 10 ) may be due 
to the way data were edited before being analyzed. Calibration corrections may 
have also been different. 

The C^(t) term in figure 11 is a measure of hemispherical or pole-to-pole 
differences. Its annual cycle is seen to have a nearly perfect sine shape, with 
a total range between its minimum and maximum values of about 20 W-m'“^. It does 
not oscillate about zero but has a bias of 2.5 W-m"^. If the Earth is symmetric 

about the equator , one would expect the terms for odd n values to have 

time histories which are symmetric about the time axis. Thus, the 2.5 W-m~2 

bias in the annual cycle of C?(t) is presumed to be due to land and ocean dis- 
tributional differences between the Northern and Southern Hemispheres. The 

C 2 (t) term in figure 11 may be considered to be a measure of equator -to- pole 
gradient. It is seen to have an average value of -25.6 W-m*"2^ with a small 
variation. 


The C^(t) terms for n = 3 to n = 12 are shown in figure 12. The 

0 o 

C 3 term is nearly sinusoidal, with a total variation of 15 W-m ^ and a mean of 

o 0 

4 W-m^-^. As with C-j , this bias is presumed to be the result of hemispheric 
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differences of land and ocean distribution. The C 4 . term has a mean of 
approximately -7 W-m“2 with variation of 5 W-m“2. its shape is not so sinu- 
soidal as the C? or C 3 histories. The C 5 to cfo terms each have a 

significant annual sine component, and the Cg and Cg terms have a bias 
such that they do not change signs. The maximum absolute value decreases with 

0 0 

increasing n until Ci 7 and C 12 / which are small and show little discern- 
ible pattern. 


Whether the lack of pattern in the computed values of €77 and C?2 
due to the nature of the atmosphere or due to the limitations of sampling and 
analysis is unclear at present. 


Effect of Grid System 

The reason for using the grid-system averaging technique is to reduce 
the computer computational burden. This is done at the expense of smoothing 
the data. One way to assess the effect of grid systems on the present analysis, 
which is based on a 5° by 5° grid system, is to compare the degree variance 
plots based on a different size grid system. This comparison was done for the 

TO 

month of August 1975 for a 5° and a 2- grid system. The computational effort 

70 2 

is 4 times greater for the 2- grid system. The results are shown in figure 13. 

2 

Out to degree 7, the degree variance plots for the two grid systems are the 
same. From degree 7 to degree 12, where their contribution is small, the degree 
variances differ only by a small amount. Thus, based on comparing degree vari- 
ances, no significant improvement is gained in going from a 5° by 5° grid system 
1 ° 1 ° 

to a 2 - by 2 - grid system. 

2 2 


CONCLUSIONS AND RECOMMENDATIONS 

Based upon the results of the analysis of the first year of data from the 
Nimbus 6 WPOV radiometer, the following conclusions regarding Earth emitted 
radiation are drawn: 

1. Degree variance plots for 12 months of LW radiation data show that the 
limit for a spherical harmonic representation of the Nimbus 6 WFOV LW data is 
12th degree. The degree variance plots also reveal that most of the power is 
in the lower degree terms. The axisymraetric (zonal) terms dominate, with their 
coefficients representing approximately 80 percent of the degree variance for 
a 12 -month period. 
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(2) Contour maps of the radiation field are computed. These maps show 
the geographical distribution of Earth emitted radiant exitance for 12 months. 
The maps reveal distinct highs and lows and their shift with the seasons. 

Highs are generally located between -35^ latitude and +35^ latitude. A band 
of lows is located between -10® latitude and +20® latitude. Lows extend to 
high latitudes in both hemispheres, with the extension being a dominant feature 
in the Northern Hemisphere. 

(3) Some of the spherical harmonic coefficients have physical meanings; 
for example, Co(t) is the global average emitted radiation, C^(t) repre- 
sents the pole-to-pole gradient, and C 2 (t) is a measure of the equator-to-pole 

gradient. Comparison of Co(t) from the present study with the global average 
of emitted radiation computed by Ellis et al. (J. Geophys. Res., vol. 83, 

0 — 

no. C4, Apr. 20, 1978) shows that CQ(t) is about 4 W-m ^ lower than their 
global average, although the annual trend is the same. The global average from 
Jacobowitz et al. ( J. Atmos. Sci., vol. 36, no. 3, Mar. 1979) shows better 

agreement with CQ(t). The C^(t) term has a bias of about 2.5 W-m“^^ pre- 
sumably due to distributional differences between land and ocean in the Northern 
and Southern Hemispheres . 

4. Based on comparing degree variances, no significant improvement is 

1 ® 1 ® 

gained in going from a 5® by 5® grid system to a 2- by 2- grid system. 

2 2 

Furthermore, the averaging of data over 5® by 5® regions is a reasonable 
approach. 

Subsequent years of data from the Nimbus 6 ERB are needed in order to 
establish a longer term average of the annual cycle and to provide a basis for 
studies of interannual variability. Nonaxisymmetr ic terms need to be analyzed 
for annual variation- Longitudinal spectra of emitted radiation as a function 
of latitude need to be computed in order to gain insight into processes causing 
a nonaxisymmetr ic distribution. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 6, 1980 
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• APPENDIX 


INTEGRATION OF ASSOCIATED LEGENDRE POLYNOMIALS 


The problem is to determine recursive formulas for Ijlj, where 


I™(© 1 ,© 2 ) = ) P™(cos ©) sin © dB 


(Al) 


where P™(cos ©) is an associated Legendre polynomial of degree n and 

order m and © is colatitude. Recursive formulas for I™ are given by 
Fricke (1978). The purpose of this appendix is to expand the derivation and 
to express the results in a form compatible with this paper. 

Equations (A2) and (A3) are obtained from Korn and Korn (1961) as follows 
(eq. (A2) is equivalent to eq. (21.8-55) in the reference): 


(n - m + 1) Pn+i (cos 0) = (2n + 1) cos © p™(cos ©) - (n + m) Pn_-| (cos ©) 


(0 ^ m ^ n - 1) (A2) 


0 


d 


\ 


[l 


p1!(COS 


>1 


(0 < m < n) (A3) 


Integration of equation (A2) using equation (Al) gives 


In+1 (®1 '® 2 '> 


(2n + 1) (n - m + 1)-'' 



rn , . 

P^(cos ©) 


cos © sin 0 <39 


- (n + m)(n - m + 1) '' ln-i(®i,©2) 


(A4) 


20 



APPENDIX 


Integration of equation (A3) gives 
» ©2 


r sin^ © 

d(cos 0) 


I^P^Ccos 0)J sin© (30 = + pJJ(cos 0) sin © cos © d© 

©1 


- (n - m 


p©2 

+ 1) ) 


Pn+-| (cos ©) sin © d5) (A5) 


Integration of the left-hand side of equation (A5) by parts gives 


r 

(n - 1) \ pJJ(cos 0) cos © sin 0 d© = -sin^ © pJIJ(cos ©) 
^©1 


©7 


©1 


+ (n - m + 1 ) l”+i (©1 ,© 2 ) (A6) 


Elimination of the integral between equations (A4) and (A6) and collecting 
terms gives the desired recursion equation 


In+1 (©1 ,©2) = 

(n + 2)*"'' (n - m + !)“'■ 

(2n + 1) sin2 0 Pn(cos 0) 

©2 

+ (n - 1 ) (n + m) Iq„i (©] ,© 2 ) 




0^ 


/m = 0, 1, 2, . . . \ 

yn = in + 1, m + 2, . . ./ (A7) 

The recursion equation (A7) for In+1 requires two initial values, 

namely, and Ini+i • Equation (A8) is obtained from Korn and Korn (1961) 

as follows (eq. (A8) is equivalent to eq. (21.8-54) in the reference): 

pJ[{(cos 0) = (1 ) (3) (5) . . . {2m - 1 ) sinJ" 0 (m = 1 , 2, , . . ) (A8) 

and from the definition of iJJ (eq. (A1 ) ) , 





(1)(3)(5) ... (2m - 1) sinn>+l 0 d0 


(A9) 
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By integration. 


= O) (3) (5) ... (2m - 1) 


-(m + 1)“^ sin"* 0 cos © 

02 

+ m(m + 1 ) “ • 

r®2 n 

\ sin"^l 0 d0 


©1 



Substituting from equations (A8) and (A9) gives 


Im(®lf®2) = -(m + l)**^ cos 0 Pm(cos 0) 


0^ 


©n 


+ m(2m - 1 ) (2m - 3) (m + 1 ) "I (©i ,0 


2 ) 


(ni — 2, 3/ ■ • .) 


Evaluating equation (A3) at n = m gives 


sin^ © fPm(cos 0)1 = (m + 1 ) cos © p!I!(cos ©) - p!II. n (cos ©) 

d(cos ©) L J ™ ">+' 


Differentiating 


m 


of equation (A8) and simplifying yields the result 


^m+1 ' (2m + 1) cos 0 pJU(cos 0) 


(m = 0, 1 , . . . ) 


Thus, equation (AlO) can be simplified as 


lIS(01,©2) = -(2m + I)-! (m + !)-■' Pm+Kcos ©) 


+ m(2ra - l)(2m - 3) (m + I)"'’ l5ri(®l'®2) 


(in — 2, 3/ • • .) 


To start this recursion, the values Iq and are needed, or 


X 02 

(1) 

©1 


sin 0 d0 = “COS 0 


©2 

©1 


(AlO) 


(All) 


(A12) 


(A13) 


(A14) 
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1 f®2 1 

I (01 ,© 2 ) = I (sin 0) sin 0 <30 = -((.■!) - sin © cos ©) 
' ‘^©1 


©2 

©1 


Next, an expression for Im+l is <3erived. From equations (A12) an<3 

r ©2 

p“+l (cos 0) sin © (33 = (1) (3) (5) . . . (2m - 1 ) (2m + 1 ) ) sin^+l 0 cos 0 
-01 '^01 


or 


Ijn+l(©l,© 2 ) = (1)(3)(5) ... (2m- 1 ) (2m + 1 ) (m + 2) -I sini”+2 © 


©2 

©1 


Simplifying with equation (A8) gives 


im+l(®l'®2) = (m + 2) ^ (2m + 3) "I pJ[|+ 2 (cos ©) 




0-j 


(m = 0, 1 , 2 , . • . ) 


m 


Thus, in suinmary, the five recursive formulas for are 


Iq (© 1 f©2) = -cos 0 


©2 

©1 


1 1 

(0-| f©2) = - (0 “ sin © cos ©) 


©2 

©1 


Im(©l*®2) = -(2m + 1)"'' (m + 1) '' p{S+i (cos 0) 


02 

©1 


+ m(2m - 1) (2m - 3) (m + I)”'' iJJJli (®1 '©2) 


(m— 2, 3f • ■ •) 


I 


(AT 5) 

(A8) , 
d© 


(A16) 

(A14) 

(A15) 

(A13) 
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(m + 2)-l (2m + 3)-1 


p;i;+2(cos 0) 


©2 

01 


(m = 0, 1 , 2, 


(AT 6) 


In+1<®!'®2) = (n + 2)"^ (n - hi + 1)"1 (2n + 1) sin^ © p”(cos 0) 


02 

+ (n - 1) (n + m) In-i (®i ,®2) 

01 


( in= 0 f 2 , • • ■ 

n = m + 1 , m + 2, . 


(A7) 
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TABLE I.- EIGENVALUES OF MEASUREMENT OPERATOR 
[tg = 6408.165 km; h = 1070 km] 



3 


h 











TABLE III,- SPHERICAL HARMONIC (NONAX I SYMMETRIC) COEFFICIENTS OF LONGWAVE RADIANT EXITANCE (W-ITl 


"cosine terms are below stairstep and sine terms are above; coefficients 
I are for 12th degree expansion at top of atmosphere; results are for limb- 
ic darkening model 


(a) July 1975 


2.996 

-.076 

-1.129 

4.639 

4.104 

-1.109 

1.706 

2.954 

-.971 -. 

-.554 

-1.382 

-2.284 -2. 

■5.142 

-3.437 

.325 -2, 

•2.326 

-1.809 

2.236 

1.706 

.367 

.403 

1.199 

1.988 

-1.799 1, 

.451 

2.770 

.034 

.318 

1.067 

1.373 

-.277 

-2.851 

.091 

. 194 

-1.889 

-1.180 


3 - 

1.0^7 -.630 .021 .263 1 . 1^9 -.103 .566 

-.045 -.951 .516 -.006 1.620 1.066 1.576 

-.222 .651 1.267 -.398 -.146 1.114 .342 


1.105 

.012 

-1.769 

-.931 

-.420 

-.215 

-1.372 

-1.946 

-.817 

.121 

-1.136 

.756 

.746 

-1.062 

1.058 

2.731 

1. C 16 

2.431 

.472 

-.588 

-.061 

-.045 

-.422 

-1.552 

.905 

-1.499 

.411 

-.923 

.508 

-1.480 

-.268 

-.472 

.298 

-.324 

.204 

-.525 

7 

8 

9 

10 










TABLE III 


(b) Augi 


1 

2 

3 

4 

5 

6 

7 

8 
9 
0 
1 


12 

11 

10 

9 

8 

7 

-.103 

-.158 

-.202 

. 327 

-. 233 

.375 

3.173 

.837 

-.351 

.329 

-.072 

-.527 

3 . 586 

4.425 

.043 

-.062 

.673 

-.187 

1.080 

2.713 

1 .124 

-.524 

.835 

.487 

- 1.134 ■ 

- 2.459 

- 1.078 

- 2. 795 

.156 

.151 

- 3.404 ■ 

- 2.701 

.230 

- 1.218 

- 1.664 

.292 

-.144 

-.075 

2.057 

. 529 

.960 

1 . C 50 

1.982 

.590 

.853 

.884 

. 749 

.321 


-.345 1.713 - 1.306 .618 

-.855 1 . C 99 - 1.443 -.245 

.116 .637 .841 -.255 

. 494 - 1.913 1 .054 .323 













TABLE III.- Continued 


(d) October 1975 



















TABLE III 


(g) Jam 




























TABLE III.- Continued 
(h) February 1976 


210 - 1.183 


1 1 . 669 | .840 -.576 -.090 .145 1 

2 .690 - 2.107 .692 -.084 -.906 

3 -319 2.229 1.780 .277 -.982 - 

4 .479 2.923 - 2.418 2.732 -.512 

1 

5 -.978 .533 .440 - 1.105 - 1.238 


6 

.5 82 

1.512 

4.196 

-.944 

7 

1 

1.748 

.094 

-.395 

1.284 

8 

.359 

- 1.253 

- 2.143 

-.524 

9 

- 1.619 

.409 

-.021 

-.139 

10 

1.116 

1.237 

-.552 

1.399 

11 

.094 

- 1.104 

-.168 

.373 


417 -.828 

797 .317 

389 .988 

628 -.498 

165 - 1.200 
?63 1 . 431 

555 1.243 

390 - 3.316 


728 

1.245 

- 1.060 

-.518 

12 

522 

1.082 

-.538 

.489 

11 

448 

- 1.366 

1.811 

2.314 

10. 

846 

-.471 

.575 

- 1.376 

9 

761 

2.456 

- 1.759 

- 1.082 

8 

062 

-.939 

-.417 

.296 

7 

201 

- 2.612 

4.009 

2 . 647 ; 

6 


270 - 


U 316 .152 1.137 1.892 -.130 5 

.010 - 2.217 1.543 - 3.014 .198 4 

.166 . 756 - 1.195 -.374 -.754 3 

205 -.442 .9191 2.976 


438 -.062 - 1.382 1.287 -.969 
725 .340 -.378 -.212 - 1.808 


/k 1 






TABLE III 


(i) Mart 



Continued 


:h 1976 


6 

5 

4 

3 

2 

1 

.761 

“.276 

“.243 

.174 

“1.539 

• 263 

“.511 

.354 

“.496 

.492 

.872 

“1.153 

.036 

.337 

.383 

“1.160 

2.918 

2.213 

051 

“1.013 

.825 

.078 

.170 

“1.136 

.385 

“1.638 

“1.652 

1.851 

“2.549 

“1.146 

.5^7 

2. 115 

“.934 

“1.295 

.077 

“1.241 

“1 .089 

1.513 

2.201 

“2.184 

4.603 

3.380 

-.80^ 

-3.260 

.041 

1.466 

1.071 

.932 

-.131 

.717 

“3.223 

.223 

“3.054 

1.096 

.5 60 

. 199 

1.122 

“1.324 

.025 

“1.316 

“.074 

“.779 

.633 

.411 

2.261 

“.481 

-.055 

.700 

.579 

“.757 

.029 

.332 

“.169 

“.506 

1.066 

.170 

“.304 

.753 
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Figure 1.- Time periods for longwave radiation measurements from Nimbus 6 (dark bands) 
for 12-month period from July 1975 through June 1976. 
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Figure 2.- Range of Sun zenith angle for which wide-field-of-view 
(WFOV) sensor measurements were deleted (111.5° to 123.5°). 
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Figure 3.- Zonal sampling plot for September, showing zonal deficiency 
in number of measurements due to removal of Sun contaminated 
measurements. 



Radiant exitance 



Orbit 


CJ 


Figure 4,- Composite plot of raw longwave radiant exitance showing some anomalies in data. 
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(b) August 1975. 

Geographical distribution of Earth emitted radiant exitance at 
top of atmosphere (W~m“2) for 12-month period. 
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(h) February 1976. 
Figure 8.- Continued. 
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(k) May 1976. 
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Figure 8.- Concluded. 
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Figure 10.- Time histories of global mean from Nimbus 6 wide-angle radiometer for July 1975 
through June 1976 (29-month composite data set for 1964 through 1971 represented by 
dashed line) . 
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